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Abstract We revisit the one-dimensional Burgers equation in the inviscid limit for white-noise initial 
velocity. We derive the probability distributions of velocity and Lagrangian increments, measured on 
intervals of any length x. This also gives the velocity structure functions. Next, for the case where the 
initial density is uniform, we obtain the distribution of the density, over any scale x, and we derive the 
density two-point correlation and power spectrum. Finally, we consider the Lagrangian displacement 
field and we derive the distribution of increments of the Lagrangian map. We check that this gives 
back the well-known mass function of shocks. For all distributions we describe the limiting scaling 
functions that are obtained in the large-scale and small-scale limits. We also discuss how these results 
generalize to other initial conditions, or to higher dimensions, and make the connection with a heuristic 
multifractal formalism. We note that the formation of point-like masses generically leads to a universal 
small-scale scaling for the density distribution, which is known as the "stable-clustering ansatz" in the 
cosmological context (where the Burgers dynamics is also known as the "adhesion model"). 

Keywords Inviscid Burgers equation • Turbulence • Cosmology: large-scale structure of the universe 



1 Introduction 

The Burgers equation which describes the advection of a velocity field by itself, with a non-zero 
viscosity, is a very popular nonlinear evolution equation that appears in many physical problems, see 
Q for a recent review. It was first introduced as a simplified model of fluid turbulence, as it shares the 
same hydrodynamical (advective) nonlinearity and several conservation laws with the Navier-Stokes 
equation. Even though it was shown later on by [2(| and [|[ that it can be explicitly integrated 
and lacks the chaotic character associated with actual turbulence, it still retains much interest for 
hydrodynamical studies, particularly as a useful benchmark for approximation schemes On the 
other hand, it has appeared in many other physical situations, such as the propagation of nonlinear 
acoustic waves in non-dispersive media (l~7j , the study of disordered systems and pinned manifolds [Tlj , 
or the formation of large-scale structures in cosmology [THHH]. There, in the limit of vanishing viscosity, 
it is known as the "adhesion model" and it provides a good description of the large-scale filamentary 
structure of the cosmic web [Hj]. In this context, one is interested in the statistical properties of the 
dynamics, starting with random Gaussian initial conditions [2ll[T9| (i.e. "decaying Burgers turbulence" 
in the hydrodynamical context). Moreover, in addition to the velocity field, one is also interested in 
the properties of the density field generated by this dynamics, starting with an initial uniform density. 

This problem has led to many studies in the inviscid limit, focusing on power-law initial spectra 
(fractional Brownian motion), Eo(k) oc k n , especially for the two peculiar cases of white-noise initial 
velocity (n = 0) @, HH [H, 0]| or Brownian motion initial velocity (n = -2) [H, [H, 1, H3 . The initial 



Institut de Physique Theorique, CEA Saclay, 91191 Gif-sur-Yvette, France 
E-mail: patrick.valageas@cea.fr 



2 



velocity fluctuations are dominated by short wavelengths in the former case and by large wavelengths 
in the latter case. Therefore, they provide two simple examples for two more general classes of random 
initial conditions, associated with —1 < n < 1 and — 3 < n < — 1 [13, [HI], which show both common 
and different significant behaviors. For instance, the integral scale of turbulence, L(t), and the tail of 
the shock mass function, scale with n as L(t) ~ t 2 /(™+ 3 ) and ln[n(> m)] ~ —m n+3 over the whole 
range — 3 < n < 1 [H, [H, El, HH, even though shocks are dense for — 3 < n < — 1 but isolated for 
— 1 < n < 1 [HI. Then, the specific advantage of these two cases, n — and n = —2, is that in both 
cases the initial velocity field is built from a white-noise stochastic field (either directly or through 
one integration), which gives rise to Markovian processes and allows to derive many explicit analytical 
results. 

In parallel with a study of the Brownian case in [13], we revisit in this article the white- noise 
case, taking advantage of the results obtained in [I3 |. In particular, we pay attention to issues that 
arise in the hydrodynamical context (velocity structure functions, Lagrangian displacement field) as 
well as in the cosmological context (statistics of the density field). Thus, the main goal of this article 
is to provide explicit results for the distributions of velocity increments and density fluctuations. As 
explained above, this complements the study [H[ of the Brownian case, so that we now have explicit 
exact results for these quantities for the two representative cases n — and n = — 2. This should prove 
useful to check the validity of approximation schemes devised for generic initial conditions and higher 
dimensions, as in 33] where the tails of these probability distributions are studied in the general case. 

We first describe in section [2] the white-noise initial conditions and the standard geometrical in- 
terpretation in terms of parabolas of the Hopf-Cole solution of the dynamics 0- Then, we recall in 
section |3] the Eulerian one-point and two-point distributions, p x (q) and p Xl ,x 2 (9i> Q2), associated with 
the inverse Lagrangian map x 1— > q, that were obtained in [131 ]. This allows us to derive in section H 
the distributions of the inverse Lagrangian increment and velocity increment, as well as the velocity 
structure functions. We also describe the limiting large-scale and small-scale distributions. Next, we 
consider in section[5]the distribution of the density within intervals of size x, and the density two-point 
correlation and power spectrum. Then, turning to a Lagrangian point of view, we study the Lagrangian 
displacement field in section [6] Finally, we describe in section [7] how the small-scale scalings shown 
by these exact results can be generalized to other initial conditions and higher dimensions within a 
heuristic approach. 



2 Initial conditions and geometrical solution 

2.1 Equation of motion 

We consider in this article the one-dimensional Burgers equation for the velocity field v(x,t) in the 
limit of zero viscosity, 

dv dv d 2 v . , . . . 

d-t +V dx=^ W1& W 

Let us recall here that in the cosmological context, the time t in the Burgers equation |T]) actually 
stands for the linear growing mode D + (t) of the density fluctuations, the spatial coordinate a; is a 
comoving coordinate (that follows the uniform Hubble expansion) and, up to a time-dependent factor, 
the velocity v is the peculiar velocity (where the Hubble expansion has been subtracted), see [Lsl. [35|. 
In these coordinates, the evolution of the density field is still given by the continuity equation (|37[) 
below, where the density p is the comoving density. If we take v = 0, that is we remove the right-hand 
side in Eq. (JTJ) , this is the well-known Zeldovich approximation [53, HH , where particles always keep 
their initial velocity and merely follow straight trajectories. The diffusive term of (JTJ) is then added as 
a phenomenological device to prevent particles from escaping to infinity after crossing each other and 
to mimic the gravitational trapping of particles within the potential wells formed by the overdensities 
18]. Of course, this cannot describe the inner structure of collapsed ob ject s (e.g., galaxies) but it 
provides a good description of the large-scale structure of the cosmic web (23[ . 

As is well known [H,Q, introducing the velocity potential ip(x,t) and making the change of variable 
ip(x,t) — —2i>hi8(x,t) transforms the nonlinear Burgers equation into the linear heat equation. This 



3 



gives the explicit solution 



v(x, t) 



dx 



with ip(x,t) — —1v\vl 



dq 



exp 



{x-qf i) (q) 



Avt 



(2) 



where we introduced the initial condition i/'o( ( ?) — ip(<lit = 0). Then, in the limit v 
descent method gives 



2z/ 

-> + the steepest 



ip(x, t) — min 

9 



^o(g) 



2t 



and w(a;,t) 



where we introduced the Lagrangian coordinate q(x,t) defined by 

(x - q) 2 



ipo(q) 



2t 



is minimum at the point q = q(x, t). 



(3) 



(4) 



The Eulerian locations x where there are two solutions, g_ < q+, to the minimization problem 
correspond to shocks (and all the matter initially between q_ and q + is gathered at a;). The application 
g h- > t) is usually called the Lagrangian map, and x i— > i) the inverse Lagrangian map (which 
is discontinuous at shock locations) 0] . For the case of white-noise initial velocity that we consider in 
this paper, it is known that there is only a finite number of shocks per unit length [28l |3||. 



2.2 Initial conditions 

In this article, we consider a white-noise initial velocity field vo(q), normalized by 

(v (q))=0, (v (qi)v (q 2 )) = D 5{q x - q 2 ), (5) 

where (..) is the average over all realizations of the initial velocity field. The velocity potential is defined 
up to a constant, and we may choose to normalize the initial potential ipo(q) by ipo(0) = 0, whence 

Mo)= f d<z'uo(g'), (M9)) = 0, (MQi)Mta)) = D <h> for 0< qi <q 2 . (6) 
Jo 

Thus, the initial velocity potential is a bilateral Brownian motion that starts from the origin. Then, 
thanks to the scale invariance of the Brownian motion, the scaled initial potential ipo(Xq) has the same 
probability distribution as A 1 / 2 i/ , o(?) I for any A > 0. Hence, using the explicit solution ([3]) we obtain 
the scaling laws 

VM) 1 ^ ^ 1/3 <K*A 2/3 ,l), v{x,t) X ^t- 1 ' s v{x/t 2 /\l), q(x,t)^ t 2 '\(x/t 2 l\l), (7) 

where = means that both sides have the same probability distribution. Thus, any equal-time statistics 
at a given time t > can be expressed in terms of the same quantity at the time t — 1 through 
appropriate rescalings. In this article we only investigate equal-time statistics, so that t can be seen 
as a mere parameter in the explicit solution Then, it is convenient to introduce the dimensionless 
coordinates, 

Q = «, X = X -, V=-, *=% C=% with 1 = {2Dt 2 )^, (8) 

*y 'V *y fy* /y* 

which express the scaling laws (J7|) (here c is the parabola height that will be introduced below in 
Eq.@). Thus, probability distributions written in terms of these variables no longer depend on time, 
and the scale X = 1 is the characteristic length of the system, at any time. On large quasi-linear scales, 
X ^> 1, density fluctuations are small and the distributions are strongly peaked around their mean, 
with tails that are directly governed by the initial conditions (but shocks cannot be neglected). On 
small nonlinear scales, X 1, density fluctuations are large (e.g., most Eulerian intervals are empty) 
and probability distributions show broad power-law regions. These behaviors will be clearly seen in 
the following sections. 



4 



2.3 Geometrical interpretation 

As is well known Q, the minimization problem ((4]) has a nice geometrical solution. Indeed, let us 
consider the downward parabola V x>c {q) centered at x and of maximum c, i.e. of vertex (x, c), of 
equation 

V x AQ) = - iS ^ + c. (9) 

Then, starting from below with a large negative value of c, such that the parabola is everywhere well 

below ipo(q) (this is possible thanks to the scaling tpo(Xq) = X 1 / 2 ip (q) which shows that ipo(q) only 
grows as Iq] 1 ^ 2 at large \q\), we increase c until the two curves touch one another. Then, the abscissa 
of the point of contact is the Lagrangian coordinate q(x,t) and the potential is given by ip(x,t) = c. 
In order to use this geometrical construction, it will be more convenient in the following to normalize 
the potential ipo by ipo(q-) = 0, where we first restrict the system to the finite interval [<7_, q + ], and to 
eventually take the limits q± — > ±oo [13fl . instead of normalizing at the origin q — as in ©. Indeed, 
this avoids making the point q = artificially play a special role. With this choice, the initial potential 
ipo {q) is a single Brownian motion that starts from the left boundary q_ . 

For the white-noise initial conditions ((5|), the process q i-^ ipo is Markovian. Then, following the 
approach of [i~3| , from the geometrical construction ([9]) one can see that a key quantity is the conditional 
probability density K XyC (cii,ipi;q2,ip2) for the Markov process ipo{q), starting from ipi at qi, to end at 
ip2 at qi > 9i, while staying above the parabolic barrier, ipo(q) > V x , c (q), for q\ < q < qi- This kernel 
was obtained in [131 ] and we recall its expression in Appendix [X] with our notations. We also derive the 
closely related kernel E xc (q 1 ,t(j 1 ;q2,ip2',q), defined in Eg. (fM?)) , which only counts among these initial 
conditions the ones that have a last excursion below V XtC +dc in the range [q, q + dq]. 



3 Known Eulerian distributions 

We briefly recall in this section the expressions of the one-point distributions, p x (q) and p x (v), of 
the Lagrangian coordinate q(x, t) and velocity v(x, t) at the Eulerian point x. We also consider the 
two-point distributions p Xl , X2 (qi, 92) and p xi ,x 2 (li: 92)- These results were already obtained in [l3[, but 
they are the basis of our computation in the following sections of the distributions of Lagrangian and 
velocity increments, from which we obtain the distribution of the matter density, and of the distribution 
of the Eulerian increment. We give more details and explicit expressions in Appendix [Bl 



3.1 One-point Eulerian distributions p x (q) and p x (v) 

To any Eulerian point x we can associate the Lagrangian coordinate q{x,t) defined as the location of 
the minimum in Eq.Q, except at shock locations where there are two (or more) contact points between 
the initial potential ipo(q) and the first-contact parabola P XlC . Since shocks are in finite number per unit 
length [3|, [28[, Eulerian points have a well-defined Lagrangian coordinate q(x,t) with probability one. 
However, note that the Eulerian position x is usually not "occupied" by the infinitesimal mass that 
was initially located at q, as all the matter is collected within shocks (thus a given Eulerian point has 
almost surely a zero matter density) (35[. Nevertheless, through Eq.Q one can derive the properties 
of the velocity field from the Lagrangian coordinate q{x,t). 

The one-point distribution, p x (q), of the Lagrangian coordinate q at point x, can be readily obtained 
from the kernel E XyC given in Eq. (|101ll . or the kernel K XyC , as shown in [l3j]. For instance, from the 
definition of E x c we can write 

Px{q)= lim / dcd0+ E xc (q^,Q;q + ,i() + ;q), (10) 

where we normalized the initial potential by ?/>o(<7-) = and we let q± — * ±00 as the size of the 
system goes to infinity, as discussed below (|9]). Thus, in Eq. ([T0| we count all initial conditions ipo(q) 
that have a first-contact point of abscissa q with a parabola V x ^ and we integrate over all possible 
heights c. We recall in Appendix IB. II the explicit expressions of Px{Q) and Px{V), in terms of the 
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scaling variables (JSj> , see Eqs. (1102j) - (|103|) and [13] . Both distributions are related through the change of 
variable X = Q + V, that expresses the second equation ([3]). Thanks to the homogeneity and isotropy 
of the system, the distribution Px{Q) only depends on the distance \Q — X\, whereas P(V) is even 
and no longer depends on X. The asymptotic behavior of the distribution of the velocity V (and of 
the Lagrangian coordinate Q = X — V), 

\ V \ » 1 : P(V) ~ . e — ,1^1-1^,-/3 (n) 

11 v ' Ai'(-wi) v ' 

shows that P(V) decreases faster than a Gaussian at large V [2, 13]. Contrary to cases where the initial 
velocity field has no ultraviolet divergence (i.e. the initial variance cr^ (0) = (vq) is finite, as for the 
case of Brownian initial velocity [3J]), the large-?; tail cannot be directly understood from the statistics 
of rare local peaks in the initial velocity field, Here, as we have recalled above, at any time t > all 
the matter has collapsed within a finite number of shocks per unit length, which merge in the course 
of time to build increasingly massive shocks within larger voids [HI, 0, [lj| • Then, the typical velocities 
observed in the system are governed by this merging process, rather than by the initial velocities of 
regular points that would not have collided yet. Nevertheless, the cubic exponential tail (fTTj) can be 
understood as follows. A structure with a large velocity v has traveled by time t over a distance of 
order x ~ vt. On the other hand, the mean velocity Vo(q) of the mass that was initially located in 
the Lagrangian interval [91,(72], of size 9 = q 2 - qi, is v a (q) = dq' v a (q')/q = (ip 2 - i>i)/l- It is 
Gaussian with a variance o$ (q) — D/q from Eq.([6]). Since momentum is conserved by the inviscid 
Burgers dynamics Q , so that the momentum of a shock is equal to the sum of the initial momenta of 
all the particles it contains, we can associate to the velocity v and the distance x = vt the Gaussian 
ht ~ e~ v2/ < {vt) ~ e~ v3t / D , where we did not write numerical factors in the exponential. This 
gives back the cubic exponential tail (fTTj) . Even though we followed shocks in the previous argument, 
in spite of the fact that they occupy a set of zero measure in Eulerian space, this still sets the tail of 
the Eulerian velocity field as the velocity of a shock located at position x is related to the local velocity 
field as ii shock = (v(x~) + v(x + ))/2, and v(x,t) has a constant slope of 1/f in-between shocks, see 
Cubic exponential tails such as (fTTI) are characteristic of probability distributions obtained for these 
white-noise initial conditions 0, P 



3.2 Two-point Eulerian distributions p Xl ,x 2 (qi, 92) and p Xl ,x 2 (vi,V2) 

We now consider the two-point distribution, p xl ,x 2 (Qii 92), of the Lagrangian coordinates {91,92} as- 
sociated with the Eulerian locations {xi,x 2 }. We take X\ < X2, which implies that q\ < q 2 since 
particles do not cross each other and therefore remain well-ordered. One needs to consider the two 
cases, i) q\ ^92, and ii) 91 =92- The first case, associated with different first-contact points, gives the 
contribution [13J 

Px u x 2 (Qi,Q2) = 0{Q2 - Qi) J(Qi - Ai) J(X 2 - Q % ) HxlxAQiM, (12) 

where 9(Q 2 — Qi) is the Heaviside function and we introduced the functions J and TL given by 
Eqs. (|103p and (|106|) , whereas the second case, associated with a common first-contact point, gives the 
contribution 

Px u x*(QiM = KQ2 - Qi) J(Qi - Ai) j(x 2 - q 2 ) e -(Q 1 -^) 3 /3+(Q 2 -x 2 )3/ 3 _ (13) 

One can check that the function TLx lt x 2 {Qii Q2), whence the contribution Pj^ x ^ (Qi, Q2), and the con- 
tribution Pj£ x 2 {Q\, Q 2 ) are invariant with respect to uniform translations of and Qi, in agreement 
with the statistical homogeneity of the system. Then, the full distribution Px 1 ,x 2 (Qii Q2) is given by 
the sum of both contributions (Tl2|) and (fT3| . Next, the two-point velocity distribution, Px lt X 2 {Vi, V2), 
is obtained from Eqs. (fl"2"]) . (fT3|) . by using Vi = Xi — Qi. 

We can note here that, thanks to the Markovian character of the process 9 1— > •^0(9), the n-point 
distributions of the velocities Vi, and of the Lagrangian coordinates qi, factorize as [131 ] 



p Xl ,..,x n (vi, -;V n ) = p Xl {vi) p(x 2 ,v 2 \xi,v 1 ) 



(14) 
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and 

(qi, ?n) = Px 1 {qi)p(x 2 ,q 2 \x 1 ,q 1 )..p{x n , q n \x n -i,q n -i), (15) 

with the transition kernels 

/ i \ Px 1 ,x 2 (vi,v 2 ) j / i \ Px 1 ,x 2 (qi,q2) 

p{x 2 ,v 2 \xi,vi) = — r — and p{x 2 , q 2 \xi, qi) = 7 — : — , (16) 

Pxiivt) PxAqi) 

that can be obtained from the two-point and one-point distributions derived above. Again, the kernels 
p{x 2 , v 2 \xi, vi) and p(x 2 ,q 2 \xi,qi) are invariant with respect to uniform translations of the spatial 
coordinates Xi and qi. However, contrary to the case of Brownian initial velocity @, [34j], the transition 
kernel does not only depend on the two relative distances x 2 — x\ and q 2 — q± (thus it also depends 
on the third distance q\ — x\). This means that the inverse Lagrangian map, x 1— ► q, does not have 
independent increments. 



4 Probability distributions of the Lagrangian and velocity increments 

We now consider the probability distributions, p x {q) andp x (w), of the Lagrangian increment, q = q 2 —q\, 
and of the velocity increment, v — v 2 — v\, over the Eulerian distance x = x 2 — x\. These distributions 
can be directly obtained from the two-point distributions (|12j) and (|13p. but they were not studied 
in previous works (except for the singular part (I108j) associated with voids). In particular, as noticed 
in the conclusion of [13| . the asymptotics of p x {v) at large v cannot be obtained in a straightforward 
manner from the estimations of their section 5, as the latter apply to the limit of large distance x at 
fixed v\ and v 2 . Next, we shall need the distribution p x (q) to derive the distribution of the overdensity 
at scale x in section [5] 



4.1 Lagrangian increment, q = q 2 — q\, and velocity increment, v = v 2 — v\ 

The probability distribution Px(Q) of the Lagrangian increment, Q = Q 2 — Q\, can be obtained 
by integrating the sum of the bivariate distributions (I12p and (fT5|) over the variable Q\ at fixed 
Q = Q2 - Qi- This gives 

/oo 
dQ 1 J(Q 1 )J(X-Q-Q 1 ) 9(Q)Ho,x(Qi,Qi + Q) + S(Q)e- Q '/ 3+ ^- x ^/ 3 , (17) 
-00 

where 9(Q) is the Heaviside function. The second term gives a contribution of the form P^(Q) = 
S(Q) P x , given by Kq. (|108[> in Appendix IB. 31 Note that Eulerian intervals with Q — also have a zero 
matter content so that P x is also the probability for an interval of size X to be empty (see section [S] 
below where we discuss the matter density field), in agreement with the result obtained in 13] for this 
void probability. We recall the properties of this distribution of voids in Appendix IB. 31 and Fig. [T2l 

In this article we are mostly interested in the regular part, P X (Q), associated with non-empty 
Eulerian intervals, which has not been studied in previous works. From the first term in expression 
(fT7| it reads as 



Pj£(Q) = 0{Q) 2V^Xe- x3/12 J 



+io ° dsdsids 2 e s (Q- x )+( s i+ s 2)x/2+( Sl - S2 ) 2 /(4X) 



(27ri) 3 Ai(si)Ai(s 2 )Ai(si - s)Ai(s 2 - s) 
dre Xr Ai(r + Sl )Ai(r + s 2 ). (18) 



We recall in Appendix [C] an alternative expression for the integral over r that appears in Eq. (JT? 
obtained in [l3j, which is useful to derive asymptotic behaviors. Thus, at large distances, X 3> 1, 
Eq. (fT8f yields the asymptotic behaviors 



X » 1 : ~ A .^ i)2 Q-V2 e -"i*-*'/ia forO<Q«X- 2 , (19) 

P|(Q) - Ai/( ^ i)2 |I/| 3 / 2 e -^H^I 3 /i2 for |V| »1 and Q» X~ 2 , (20) 
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Fig. 1 (Color online) Left panel: The prob abil ity distribution Px(Q) of the Lagrangian increment Q, for three 
Eulerian sizes, X = 0.5, 1 and 2, from Eq. (|18fl . We have Q > and all curves display an inverse square root 
singularity oc l/\/Q at Q — *• + . In addition, there is a Dirac contribution, PxS(Q), with the weight P% 
displayed in Fig. 1121 Right panel: Same as left panel but on a logarithmic scale. 

where V — X — Q is the dimensionless velocity as in ©. At small distances, X <C 1, Eq. fTS)) leads to 

X<1: P|(Q) ~ ^=Q~ 1/2 forO<Q<l, (21) 

P^{Q)^2^XQ r ° /2 e- UJlQ - Q3/12 for 1 « Q < A~ 1/2 , (22) 

P|(Q) ~ 27iVA 7 Q 3/2 e-" lQ - Q3/12 for(J»r 1/2 . (23) 



Thus, at all scales X the distribution Pjf(Q) displays an inverse square root tail at low Q. At large 
X this tail has an exponentially small weight, that scales as the weight PS- of the empty cells, and it 
is restricted to very low Q, whereas at small X it describes the full low-Q regime. As expected, we 
can check from Eq. (f20|) that on large scales, X>1, the Lagrangian increment is centered on X, with 
the usual cubic exponential tails encountered for this white-noise initial velocity spectrum, whereas on 
small scales, I < 1, the distribution shows a monotonous decline. 

We display in Fig.[T]the probability distribution Px(Q) for three Eulerian sizes. This clearly shows 
the change of shape as we go from large to small scales, as well as the translation of the mean (Q), 
that follows X from the conservation of matter (see Eq. ([38|) below). Note that for numerical purposes, 
in order to follow the evolution of Px(Q) with X, and its behavior over the different characteristic 
domains listed in Eqs. (fTT)|) - (l2"3"l) . it is useful to gradually move the integration contours in the complex 
plane of Eq. (fT5)) as one goes from one regime to another one (but making sure that one does not cross 
singularities). One interest of these results is to provide an explicit example that is representative of 
initial conditions in the range — 1 < n < 1, where n is the slope of the initial energy spectrum (i.e. 
Eo(k) oc k n+1 ~ D in D dimensions), which show significant power at high wavenumbers. Then, we 
can see that the distribution P^(Q) always diverges at low Q (i.e. at low density) as l/y/Q. This 
implies in particular that, contrary to the cases —3 < n < — 1, the very low-Q part of the distribution 
cannot be estimated through steepest-descent approaches that apply to rare events, as discussed in 
f33l ]. Nevertheless, these approaches can give an estimate of P^(Q) in the quasi-linear limit X — > 00 
at fixed Q, studied in section B~2l below. 

The probability distribution, Px(V), of the velocity increment, V = V2 — Vi, is obtained from the 
distribution Px(Q) by using the relation V = X — Q. We show our results in Fig. [21 for the same 
three Eulerian scales as for Px(Q) displayed in Fig.[T] There is no longer a translation of the typical 
velocity, since (V) = for any scale X, but we clearly see the translation of the inverse-square root 
tail, ~ 1/ \JX — V, that follows the upper bound V < X associated with empty cells (as seen from the 
relation X = Q + V and the constraint Q > 0). 




Fig. 2 (Color online) Left panel: The pro bab ility distribution Px(V) of the velocity increment V , for three 
Eulerian sizes, X = 0.5, 1 and 2, from Eq, (|18fl . We have V < X and all curves display an inverse square root 
singularity oc 1/y/X — V at V — * X~ . In addition, there is a Dirac contribution, Pxd(X — V), with the weight 
P x displayed in Fig. 1121 Right panel: Same as left panel but on a semi-logarithmic scale. 



4.2 Asymptotic distribution on large scales 

On large scales, X ^> 1, the distributions Px{Q) and Px(V), except for the exponentially small 
contributions associated with Q = and Q <C X~ 2 in (|108|) and (fT9| . can be described by the 
symmetric distribution f^, 

X>1: P x {Q)~Foo(X-Q) and P X {V) ~ (24) 

with a Fourier transform .Foo given by: 

" /__ S e wi,t - CL Si MWTi)) ' (25) 

Equation (|25p is obtained from Eq.(|18p by using the asymptotic behavior of the integral over r, as 
given by the first term in Eq. ljllip . and next making the change of variable s — ik. The scaling function 
Foo(V) no longer depends on X: the distribution of the velocity increment V converges to the finite 
distribution on large scales X — > oo. Note that in this limit the upper boundary on V, V < X, 
associated with the positivity of Q — X — V, goes to +oo so that the limiting function J-oc(V) is 
defined over the whole real axis. This is also why the Laplace transform (fT5)) naturally gives rise to the 
Fourier transform (f25|) in this limit. Moreover, we can see from Eq. ([25|) that Foo{V) is even (the change 
of integration variable s' ~ + s' — ifc/2 in Eg. ([25)) readily shows that T x (k) is even). The asymptotic 
behaviors of TooiV) at large V can be read from Eq. (f2"Dl) : 

^(V)- - . 7~~ |y|3/2 e -^|V|-|^| 3 /12, (26) 

Ai'(— wi)^ 

Thus, we recover the cubic exponential tails that are characteristic of white-noise initial conditions 
0,1, [13 and can be understood at a qualitative level following the discussion below Eq. fTTj) . In fact, 
as shown in [33], on a quantitative level it is possible to obtain the factor 1/12 in the exponential 
(|26p through a simple steepest-descent approach, that identifies the initial conditions (i.e. the relevant 
saddle-points) that give the main contribution to these tails. Then, the exact result ([2l)|) provides 
a useful non-trivial test of such general approaches that rest on some additional assumptions (for 
instance, one only looks for symmetric saddle-points). 

Of course, the scaling function J-oo(V) does not capture the low-Q power-law tail (|19p . However, 
using the property 

r+MO ds 1 

2 7 riAi( S )2- 1 ' (27) 
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Fig. 3 (Color online) Left panel: The asymptotic di stri bution .Foo(V) of the velocity increment V , reach ed in 
the limit of large Eulerian distance X 1, from Eq. (|25[) . The dashed lines show the asymptotic behavior (I26|) . 
Right panel: Same as left panel but on a logarithmic scale. 



we can sec that the distribution J-oo(V) is normalized to unity as it should, since the weights of the 
low-Q tail (fl9|) and of the Dirac term (|108p vanish in the limit X — » oo. We show the scaling function 
Foo{V) in Fig. El 

From Eqs. (|24|) -([25ll. the moments of the velocity increment are given in this limit by 

X^+^: (V 2n+1 )=0, (^ = (-1)^(0), (28) 



dk 2n 



whence 



2 r +ioc ds s 

(V) = and (Q) = X, (V») = (Q 2 ) c = -- / ^ * 0.837 (29) 

We can check that we recover (Q) = X, as implied by the conservation of matter. We can see that 
there is negligible power on large scales since (Q 2 ) c goes to a constant for X — > +oo. This holds for 
cumulants of all orders, as the scaling function Foo(V) does not depend on X (see also the left panel 
of Fig. [B] below, where we can see that both (Q 2 ) c and (Q 4 ) c have a finite nonzero large-scale limit 
whereas (Q 3 ) c vanishes by symmetry of J 7 QO ). 

Equations (f2~i]) - (l2l)l) and (|2"5)) show that, because of the lack of power at large scales in the initial 
velocity field, at any time t > the system observed at any scale x, whatever large, is governed 
by nonlinear effects and exhibits strongly non-Gaussian statistics, even though the initial conditions 
are Gaussian. Indeed, the redistribution of matter within a series of discrete shocks has regularized 
the initially singular white-noise velocity field, through the balance between the infinite different sign 
velocities of neighboring particles, over lengths of order (2Di 2 ) 1 / 3 , that have merged in a single shock. 
Moreover, the velocity field in the voids is governed by the motion of the boundary shocks, since from 
Eq.([3]) it has a constant slope 1/t in-between shocks, and the velocity of a shock satisfies v shock = 
(v(x~) + v{x + ))/2, see [7|. These processes are clearly non-perturbative and give rise to the non- 
Gaussian statistical properties described above in the large-scale limit. This would not be the case for 
initial conditions with significant initial power on large scales. Then, even though shocks may have 
formed as soon as t > 0, one still recovers the initial Gaussian statistics on large scales, as explicitly 
checked in [3J] for the case of a Brownian initial velocity field (where the initial energy spectrum is 
Eo(k) oc ~kr 2 instead of the constant spectrum associated with the white-noise initial condition studied 
in the present article). 

Again, these exact results provide a useful confirmation of the results obtained by approximate 
methods, such as the steepest-descent approach of [Hj]. Indeed, there it is found that for n > D — 3, 
which includes the case {n = 0, D = 1} studied in this article, the relevant saddle-points always give 
rise to shocks, which is not the case for initial conditions with less initial power at high wavenumbers 
(such as {n = —2, D — 1}, i.e. Brownian 1-D initial velocity). 



10 



4.3 Asymptotic distribution on small scales 

On small scales, apart from the Dirac contribution (1108j) . associated with empty cells, and the very 
large-Q tail (f23|) . the distributions Px{Q) and Px(V) can be described by the function Tq, 

P X (Q) ~ X F (Q) with Q>0 and F (Q) = / — e sQ ^ (s), (30) 

J-ioo 27Tl 

with 

- = /-+ 100 dV 1 9 AiV + s ) = /- +loc ciV AiV)AiV + «) 
0U 7_ loo 2m Ai(s') 2 9s' Ai(s' + s) J_ loc 2iri Ai(s') 3 Ai(s' + s) ' 1 j 

The expressions (|3U)) - ([3T|) are obtained from Eq. (fT8"]) by taking the Gaussian integration over S2 
(since |s2 — Si| ~ Vx we can set at leading order S2 = Si in (JT5J), apart from the Gaussian fac- 
tor e( Sl_S2 ) /( 4X )), and next setting X — 0, which allows to perform the integral over r. Note that a 
change of variable and an integration by parts allow to write J-q{Q) as 

MQ) = 2Q \J_^ - ^ j - e — j . (32) 

It is clear that the integral over Q of the distribution (|30"|) is not normalized to unity since it decreases 
as X at small Eulerian distance. Indeed, in this limit almost all Eulerian cells have a zero Lagrangian 
increment (whence a zero matter density), associated with the Dirac contribution (|108|) (see the first 
limit in (|109p ). whereas non-empty cells occur with a probability proportional to X, see Eqs.(|21 )) - (|23)) . 
This can be directly understood from the fact that all the matter is condensed into discrete shocks 
that occur in a finite number per unit length (28l . Q , so that the probability for an Eulerian interval 
to contain at least one shock (which is equal to 1 — P x ) scales as X for small cell size X. 

In fact, the comparison of Eq. (|32"l) with results obtained in [l3j shows that Tq(Q) is also the mass 
function of shocks, as we shall check in section 16.31 below through a different method. Therefore, the 
expression (|30D actually means that on small scales the probability distribution Px{Q) is asymptotically 
equal to the probability to encounter one shock of strength Q in the interval of size X. Indeed, since 
shocks are isolated it is clear that in the limit of small size X the probability to have two or more 
shocks within X goes to zero faster than X, so that Px{Q) is governed by the probability to encounter 
one shock over the length X, which directly gives the scaling (f30|) where !Fo(Q) would be defined as 
the shock mass function. Thus, how results explicitly show how the scaling (|30|) and the shock mass 
function J-q(Q) arise from the full distribution (jT8j) of the Lagrangian increment Q. 

This property is well-known to hold for any Burgers system without dense shocks fl5l l30l|. Then, 
as pointed out in [30] , who studied the case of compactly supported white- noise initial velocity, the 
statistics of the velocity field (whence of the Lagrangian increment Q) at scales much smaller than the 
average distance between shocks are fully determined by the one-point distribution of shock strength 
n(m) (which in our case is equal to Fq(Q) as seen in Eq. ([77|) below). As noticed in [3(|, the case of 
compact initial conditions is in a different universality class than the system studied here, where the 
white-noise initial velocity (|5|) extends to the whole real line, in the sense that the scalings ((8]) no 
longer hold. Indeed, the size L of the initially non-zero velocity field introduces a new scale and at 
late times there are only two shocks left, which allows [301 ] to compute both single- and multiple time- 
velocity structure functions. Nevertheless, on small scales, much below the typical distance between 
shocks, both systems show the same scalings (see Eqs.(|35]) and (|36|) below), governed by the one-shock 
contribution. 

Again the asymptotic behaviors can be read from Eqs. ([2T|) -([23 |) : 

Q«l: F {Q)~-^=, Q»l: F (Q) ~ 2^Q 5/2 e -^ Q - Q ^ 12 . (33) 

We show the function !F (Q) and its asymptotic tails (|3"3"|) in Fig.[H see also [HI]. Note that the scaling 
function Fq{Q) does not describe the very far tail Q ^> X^ 1 / 2 of Eq. (|2"3"|) , which is repelled to infinity 
in the limit X — >■ 0. This very high-Q tail is related to the behavior of P X (Q) at large scales, as seen 
from the comparison with Eq. (|23p . Indeed, it corresponds to the limit of very rare events, where the 
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Fig. 4 (Color online) Left panel: The scalin g fun ction J-q(Q) that describes the distribution of the Lagrangian 
incr eme nt Q in the limit I < 1, from Eqs. (|30[) - (|31[) . This is also the mass function of shocks, as checked in 
Eq, (177[l below. The dashed lines are the asymptotic behaviors (|33[) . Right panel: Same as left panel but on a 
logarithmic scale. 



tail of the distribution is governed by specific initial conditions, independently of the scale X. These 
are the saddle-points obtained in (33[, which set the cubic exponential falloffs of both Eq. and 
Eq. ([23[) . Thus, for any finite X the very far tail ([23]) of the distributions Px(Q) and Px(V) is not 
captured by the shock mass function, but this regime is repelled to infinity as X — > 0. 

For any v > 0, where the contribution from the Dirac term (| 1 08[) vanishes, we obtain for the 
moments of the Lagrangian increment 

r+ ioc ds 

u>0: (QVirKl] HO""- ^oW, (34) 

where the integration contour runs to the left of the origin, Sft(s) < 0. Thus, we recover the fractality of 
the inverse Lagrangian map, (Q v ) oc X, which is well known to be due to the contribution from shocks 
as discussed above [HI, [301 ] . Indeed, if we have a shock of finite Lagrangian length SQ S at position 
X s , it gives a contribution [Q(X S + X/2) — Q(X S — X/2)] v ~ {5Q S ) V which remains of order unity 
for X — > + for any v > 0. Next, the probability to have a shock of a given finite strength SQ S in a 
small Eulerian interval X scales as X at small distances, which gives rise to the factor X in Eg. ([34]) . 
Therefore, the scaling (|34l) is actually quite general and applies as soon as shocks have formed with 
a finite probability [la . 1301 ] , above a critical exponent v c that depends on the initial conditions (here 
v c = 0). We can note that the moments diverge for v < because of the Dirac contribution p08[) . 
whereas for other initial conditions such as a Brownian initial velocity they can remain well-defined 
and obey a second scaling law below v c [HI, [l| . For integer v we obtain: 



n > 



1: (Q»)~X(-l)"-^(0), whence (Q) ~ X, (Q 2 ) - X— J 1.136 X. 

(35) 

Again, we can check that (Q) = X, in agreement with the conservation of matter. The scaling (|35[) 
also implies for the cumulants (Q n ) c «Xin the small-scale limit X <C 1, as can be checked in the left 
panel of Fig. [S] This gives for the moments of the velocity increment 

(V)=0 andforn>2: (V n ) ~ X -^(0), ([«(x a , i)-«(x l5 1)]») ~ ^ l) (0) f— J ^"j^ - 

Thus, we recover the usual anomalous scaling of the structure functions, ([v(x + i) — v(x)] n ) oc £ at 
small distance t, that was also observed in numerical simulations [28[. As explained above, this is due 
to the contribution from shocks [3, [3(| • We further discuss these small-scale scalings in section 17.11 
below, on a more general setting. 
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?7 = m/(p x) log(r;) 

Fig. 5 (Color online) Left panel: The probability distri buti on Px fjf) of the overdensity n = m/(pox) = Q/X, 
for three Eulerian sizes, X — 0.5, 1 and 2, from Eqs. (|18[) , (|38|l . All curves display an inverse square root 
singularity oc 1/y/rj at small densities -q — > + . In addition, there is a Dirac contribution, P^5(rf), with the 
weight P x displayed in Fig. 1121 associated with empty cells. At large scales X we recover a distribution that 
is sharply peaked around the mean density, (rj) = 1 (i.e. (p) = po). Right panel: Same as left panel but on a 
logarithmic scale. 



5 Density field 

5.1 Distribution of the overdensity r\ as a function of scale x 

We now consider the evolution of the density field, p(x,t), that is generated by the Burgers velocity 
field, starting at t = with a uniform density pq. Thus, the density field obeys the usual continuity 
equation, 

^ + ^(H = 0, with p(x,0)=p , (37) 

whereas the velocity field evolves through the Burgers equation (|TJ) . As recalled in the introduction, in 
the cosmological context this also provides an approximation for the formation of large-scale structures 
(the cosmic web), known as the "adhesion model" [H,|3^|. Then, the density p is the comoving density 
(i.e. measured in comoving coordinates x that follow the Hubble expansion) and t is linear growing 
mode D+. Thanks to the conservation of matter, the mass m located between the Eulerian positions 
Xi < X2 is m = po{q2 — qi), where q(x, t) is the inverse Lagrangian map. Then, the overall overdensity, 
i] = m/(pox), in the interval of size x — X2 ~ xx, is given by i] — q/x = Q/X. Thus, the probability 
distribution, Px(t]), of the overdensity 77, is given by the distribution, Px(Q), of the Lagrangian 
increment Q, through 

ri= — = % whence P x (r}) = XP X (Q). (38) 
Pox X 

Explicit expressions are obtained by substituting Eqs. (|108l) and (|18p. In particular, the singular part 
(|108p gives the Dirac contribution P x 5(rj) associated with empty cells. For the regular part, r\ > 0, the 
asymptotic behaviors at large and small scales are directly read from Eqs. (fTi?)) - (|23[) as 

X » 1 : P x fa) ~ . , * . g q-Va e-" 1 *-* 3 / 12 for v « X~ 3 , (39) 
Ai'(— u>x) 

p x(v) ^ /f X5/ " 2 \V - 1| 3/2 e -^W-i|-^ 3 W-H 3 /i2 for 1 _ !| » x -i and » x -3 5 (40) 

and 

AT 3 / 2 

A<1: P X M ^V~ 1/2 for? ? <A- 1 , (41) 

P X ( V ) ~ 2V^A 9 / 2 ry 5 / 2 e -^n-x^ 3 /l2 for x -i <<; ^ <<; ^-3/2^ (42) 

P X ( V ) ~ 2^A 3 ry 3 / 2 e -^x^-^ 3 ,7 3 /l2 for 7] » ^-3/2 (43) 
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Fig. 6 (Color online) Left panel: The first few cumulants {Q n )c of the Lagrangian increment Q, as a function 

of the dimensionless scale X = x / '(2Dt 2 ) 1//3 . Right panel: The ratios SV,, defined by the first equality in (|48[) , 
as a function of X on a semi-logarithmic scale. 



Thus, at all scales we have an inverse square root tail at low densities, oc 1/y/rji that follows from 
the low-Q tail obtained in section E3~T1 Again, its weight shows the same cubic exponential decay at 
large scales as the weight P x of empty cells (I109p . whereas at small scales it describes the full low- 
density regime. On large scales, the density distribution is centered on the mean (77) = 1, with cubic 
exponential tails on both sides, until it reaches the very low-density tail oc 1/y/rj at 77 <C l/X 3 . On 
small scales the density distribution shows a monotonous decline, with again a cubic exponential tail 
at large densities, 77 3> l/X. We display the density distribution Px(v) m Fig. El for three Eulerian 
sizes as in Fig. [TJ to show its evolution with scale. Again, the cubic exponential tails (|40|) and (|43j) can 
be obtained from a simple and general steepest-descent method [33J. However, the exponent —1/2 of 
the power-law regime that appears at small scales in Eq. (|4Tj) is beyond the reach of such methods. It 
would be interesting to build general approaches that would be able to describe this highly nonlinear 
regime, for generic initial conditions and dimensions. Then, the result (I41|) would allow one to check 
the accuracy of such a method for a non-trivial case. 

From the results of section 14.21 we can see that on large scales the distribution of the overdensity 
is described by the asymptotic distribution 

A»l: P X (7 1 )^XT QC (X{ V -1)), (44) 

which is increasingly peaked around 77 = 1 at larger scales. Thus, we recover as expected the uniform 
density po on large scales, with a distribution that falls off faster than a Gaussian, as e~ x l r?_1 I 12 . 
Since 77 = Q/X = 1 — V/X, we obtain from Eqs. ([2"8]) - (|2"9j) for the moments of the density in this 
large- A limit: 



A » 1 : ((77 - 1) 2 " +1 ) = 0, ((77 - l) 2 ") = i-JL ^ 2n) (0), whence (77) = 1, ( V 2 ) 



OO V"/! y,, ^, y, /C — ^ 



°- 837 (45) 



In agreement with the conservation of matter we can check that (77) = 1. 
On small scales we obtain from section 14.31 



A«l: P x {r,) ~ A 2 ^ (At7), and for v > : (rj") ~ A~' y+1 r[v + 1] / ^ (-s)^ 1 ^ (s). 

J-ioo 27T1 

(46) 

As seen in section POl the scaling function Tq{Q) in (|30[) describes the probability distribution of 
the Lagrangian increment down to Q = + , hence it also describes the probability distribution of the 
overdensity down to 77 = + in (|46|) . Next, the moments of integer order are given by 

A«l, 77>1: (77") ~ A-" +1 (-l)"^ n) (0) whence (77) = !, ( V 2 ) ~ H^, (47) 
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which gives the cumulant hierarchy 
X < 1, n > 1 



Sn (X) = ~ (-1)" (Q ) , (48) 



Thus, the ratios S n (X) have a finite limit for X — > and the cumulant generating function ifx(y) can 
be written as 



X « 1 : ^(y) = ^(-l)- 1 ~ -J^'(O) 



n=l 



( -J— ) - Mo) 



(49) 



As explained in section [4~3l the scalings (|46]) - (l49p are due to the presence of shocks and are therefore 
quite general: they apply as soon as shocks have formed with a finite probability, for any initial 
conditions. We discuss these small-scale scalings in sections 17.21 and T7.3I below, for more general initial 
conditions and higher dimensions. 

We can note that on large scales the ratios S n (X) go to zero for odd n and diverge for even n, 
as seen from (l4"5"1) . This is due to the fact that, even though we start with Gaussian initial conditions 
at t = 0, the initial energy spectrum is so "blue" (which also leads to a singular white-noise initial 
velocity) that at any time t > the system is strongly affected by nonlinear effects (associated with 
the building of isolated shocks amid empty regions). This regularizes the density distribution, p x {il)i 
but the latter remains non-Gaussian in the large-scale limit x — > oo, as seen in Eqs. (pF4"|) - (H5|) or the 
explicit expression (40]). Thus, we have (r] 2n ) c - {rj 2 )™ and S 2n {X) ~ (rj 2 ) 1 c ~ n ~ X 2(n - 1 \ for n > 1 
and X — > oo, whereas odd cumulants are exponentially small. By contrast, for initial conditions with 
a sufficiently "red" spectrum, such as the Brownian case, the density distribution becomes Gaussian 
on large scales (it remains governed by the initial field and linear theory) and the ratios S n (X) have 
a finite large-scale limit that can be computed through perturbative means, see [H], H3|. 

We show in Fig.[6]the first few cumulants {Q n ) c and ratios S n - As explained above and in section l4~2l 
at large scales both (Q 2 ) c and (Q 4 ) c reach a nonzero value, as the large-scale distribution is non- 
Gaussian, whereas the odd cumulant (Q 3 ) c shows a cubic exponential decay, since the limiting scaling 
function Too is even. At small scale all cumulants show a linear dependence on A, in agreement with 
(|35| . Then, the ratio S3 also shows a cubic exponential decay on large scales whereas S4 goes to —00; 
next on small scales both coefficients reach a nonzero asymptotic value. 

5.2 Density two-point correlation and power spectrum 

We now consider the two-point correlation, £(x, t), of the density field p(x,t) itself: 

(p(xi,t)p(x 2 ,t)) c = pl£{x2 - xi,t), whence (r] 2 ) c =[ (ixi ^ x ' 2 £( X2 _ Xl y (50) 

Jo x 

In terms of the dimensionless variables ([5]), using r\ — Q / X, we obtain 

(Q 2 )-X 2 = J^ dA 1 dX 2 C(A 2 -A 1 ), whence £(X) = \^(Q 2 ) ~ 1- (51) 

Then, the small-distance behavior (|35|) . that was associated with shocks, gives rise to a Dirac contri- 
bution 

?(X)=J%(Q)S(X), whence f(x) = {2Dt 2 ) 1 ' 3 J% (0) 5{x) ~ 1.136 (2Dt 2 ) 1 / 3 S(x). (52) 
In terms of the density power spectrum, V(k, t), defined by 

/°° dk 
— e ikx p(k,t), (p(k u t) P (k 2l t)) =5(k 1 +k 2 )2TrV(k 1 ,t), (53) 
-00 27T 

this gives the asymptotic behavior at high k, 

fc->oo: V(k,t) T '(0)(2Dt 2 ) 1 / 3 ~ 1.136 {2Dt 2 ) 1 ' 3 . (54) 
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Fig. 7 (Color online) Left panel: The density two-point correlation, £(x) = {p(xi)p(xi + x)) c /po, as a function 
of the dimensionless scale X = x/(2Dt 2 ) 1 ^ 3 , fr om Eq. (|56[) . It is negative over x > but there is an additi ona l 
Dirac contribution at the origin, given by Eg . (152fl . The dashed line is the large-X asymptotic behavior (|58l) . 
Right panel: Same as left panel but on a logarithmic scale. 



As expected, shocks, that form a series of Dirac peaks in the density field, give rise to a white-noise 
power spectrum in the limit of high wavenumbers. In addition, there are also non-zero correlations 
at finite distances, that can be obtained from the second cumulant of the Lagrangian increment Q 
through (|5Tj) . Using the expression (fl8|) . making the change of variable Si — > s, ; + s and integrating 
over Q and s, we obtain 



2 n/ -^_ x3/12 [ + ^d Sl d S2 e (si+S2)X/2+(sis 2 )-/(4X) 

\Q / — 



ds 2 



2VwXe 

o 



(27ri) 2 Ai(si)Ai(s 2 )Ai(si + s)Ai(s 2 + s) 
dre Xr Ai(r + si + s)Ai(r + s 2 + s). (55) 



Next, going back to Si — ► Si — s, using the expression (|111| . which allows to remove the asymptotic X 2 
behavior of (Q 2 ), and differentiating twice with respect to X, the density two-point correlation reads 
as 

/+200 J J -f f f-OO 

-too (2tti) 2 Ai(si)Ai(s 2 ) I J x 
+ e -*.!.. a W [2< iS2 (X)/i Sl , S2 (X)+A sl , S2 (X)(^ liS2 (X)-^ liS2 (X)/ lsiiS2 (X))] }, (56) 



where the functions <? SliS2 (A) and /i SliS2 (A) are defined in Eqs. (|1 12|) - (|1 13[) and we introduced the 
function A SltS2 {X) given by 



This yields the asymptotic behaviors 



e -sX 



Ai(si - s)Ai(s 2 - s) 



(57) 



m * -0.56, and for X - co : £(A) ~ 1R .,,/ 7r A 11 / 2 e -<-^-^/ia. ( 58 ) 

16Ai'(— tJi)- 2 

Thus, as we can check in FigJTl the density correlation is negative for x > 0. This may be understood 
from the fact that, since the matter collapses within isolated zero-thickness objects (shocks), close to 
a shock there is a relative underdensity as matter has already fallen into that shock. In terms of 
particles of infinitesimal mass, the massive aggregate associated with the shock has swept matter from 
its neighborhood along its motion at previous times as particles stick together after collisions. Thus, 
starting with a white-noise initial velocity which shows no correlations over finite distance x > 0, 
some (anti-)correlations appear as soon as t > over scales of order x ~ (2_Dt 2 ) 1 / 3 (i.e. X ~ 1), but 
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Fig. 8 (Color online) Left panel: The dimensionless density power spectrum V(K), as a function of the di- 
mensionless wavenumber K = (2Dt 2 ) 1 ^ 3 k, from Eq. (|59[l . It goes to zero as K 2 at small K and it goes to a 
constant at large K. Right panel: Same as left panel but on a logarithmic scale. 



they remain very weak as they decay even faster than a Gaussian at larger scales. Again, the cubic 
exponential falloff (|58p can be understood from simple arguments. Following the discussion above, 
correlations at scale x arise from the motion of shocks over distances of order x and the building of 
voids of size X. Then, as discussed below Eq. (1109|) . this can be associated to an initial mean velocity 
vq{x) ~ x/t over the interval x and to a probabilistic weight ~ e - ' 2 ^*) I^PI X > ~ e~ x ^ Dt \ which gives 
back the cubic exponential tail (|55|) . In agreement with this discussion, we can check that the tail (|58D 
is the same (apart from the power-law prefactor) as the one obtained in Eq. (ll09|l for the probability 
of voids. 

From Eq. (|53p the density power spectrum can be written, in terms of dimensionless variables, as 



V{K) 



dX cos(KX) £_(X) + J^'(O) 



dX [cos(KX) - i]t;(x), 



(59) 



where K = jk (and 7 was defined in Eq.®). In the first equality we explicitly separated the Dirac 
contribution ([5^| from the integral over X > 0. The second equality follows by noticing that V(0) = 0. 
Indeed, from Eq. (j5Tj) we have 



V(0) 



dX£(X) = lim 



x- 



dX 



(Q 2 



2X 



0. 



(60) 



as the term in brackets decays as e - x3 / 12 at large X, as seen from Section l4~2l and Appendix [Cl We 
show in Fig.[5]our results for the power spectrum, using Eqs. ([55|) and ([59^) . We clearly see the quadratic 
behavior at low K, that can be obtained by expanding the cosine in Eq. (l59p . and the saturation at high 
wavenumbers to the white-noise spectrum (|51| due to shocks. Since the correlation function decays 
faster than a Gaussian at large distances the power spectrum is actually regular at k = 0. The high- 
wavenumber behavior (|54[) is universal and appears as soon as shocks have formed, along with the 
scalings (|35j) -(|36 |) for the Lagrangian and velocity increments observed on small scales. The quadratic 
low-fc behavior applies from Eq. (|59|) to initial conditions such that the linear power on large scales 
decays faster than k? (by contrast, if the initial velocity field is given by a Brownian motion, which 
shows significant power on large scales, the density power spectrum is exactly a white-noise spectrum 
over all k, that is "P(fc) is constant down to k = 0). 

Again, the results (|56|l and (|59|) . shown in Figs. [7] and [H may be useful to test general approxima- 
tion schemes. In particular, in the cosmological context, the matter two-point correlation and power 
spectrum are among the main observables used to constrain cosmological scenarios (both the global 
cosmological history, through the linear growth factor of density fluctuations, and the primordial initial 
conditions, generated by an hypothetical inflationary stage, through the shape of the power spectrum). 
For gravitational systems of this sort, no good approximation scheme has been obtained yet that is able 
to estimate the density power spectrum in both linear and nonlinear regimes (i.e. from large down to 
small scales), so that one needs to use numerical simulations. The case studied in this article provides 



17 



a rare hydro-dynamical example, closely related to 1-D gravitational dynamics as recalled above, where 
a complete exact solution can be derived. In this respect, the present case of white-noise initial velocity 
is somewhat more interesting than the case of Brownian initial velocity, where the power spectrum is 
simply a constant over all scales, as it shows a transition between different low- and high-wavenumber 
regimes. 

A key difference between Burgers dynamics and gravitational systems (and real turbulence) is that 
the high-fc regime is quite simple and universal, since it is governed by shocks and shows a constant 
white-noise asymptote as in Eq . ([54")) . By contrast, in 3-D gravitational (or Navier-Stokes) systems, 
small-scale structures may show a broader variety (extended halos, vortices, ...) (36l. l22l| and it is not 
known whether universal exponents exist and for which class of initial conditions they hold (in the 
cosmological context numerical simulations suggest that there is no such universality as the high-A; 
slope seems to depend on the initial slope HH). 



6 Lagrangian displacement field 

6.1 One-point distribution 

We now consider the Burgers dynamics from a Lagrangian point of view, as opposed to the Eulerian 
point of view described in the previous sections. Thus, labelling the particles by their initial position 
q at the initial time t = 0, we follow their trajectory x(q,t). Since particles do not cross each other 
they remain well-ordered. Then, the probability, p q (> x), for the particle q to be located to the right 
of the position x, is equal to the probability, p x (< q), for the Lagrangian coordinate q(x) associated 
with position x to be smaller than or equal to q. This yields 

Pq(X) = —^Px{<Q) = -^J Q dQ'P(Q'-X) = P(Q-X), whence P Q (X) = P X (Q), (61) 

where we used from Eq. (|102p the property that the Eulerian distribution Px(Q) only depends on the 
relative distance Q - X as P X (Q) = P(Q ~ X) with P(V) given by Ea.([T02"]). Thus, the one-point 
Eulerian and Lagrangian distributions are identical. This applies to any initial conditions which are 
statistically homogeneous and isotropic, so that Px{Q) only depends on \Q — X\. 



6.2 Two-point distribution and relative distance 

We now investigate the two-point distribution of the Lagrangian displacement field. In a fashion similar 
to the one-point distribution, we can relate the Eulerian and Lagrangian distributions by 

Pqi,q 2 (> X u < X 2 ) = P Xl .x 2 (< Qu> Qz). (62) 
Then, for Qi < Q 2 the Dirac part (fT3")l does not contribute and we obtain from Eq. (fT2")) . with X\ < X2, 



/Wl /-oo 
dQ[ / &Q l 2 J{Q l 1 -X 1 )J(X 2 -Q l 2 )Hx 1 ,xAQ'i,Q2)- (63) 
-OO JQ-2 

Using Eqs.dinSl) and (TiHe]) . this yields 

PQi,Q2{Xi,X 2 ) = - 2^tt(X 2 -Xi)e y 1,1 / — — -; -— ; tt- 

dXxdX-i J_ loa (2tti) 4 (si - s 2 ){s' 1 - s' 2 ) 

p -s 1 X 1 +s' 1 X 2 + (s2-s' 2 )(X 1 +X2)/2+(s2-s' 2 ) 2 /(4(X 2 -X 1 )) poo 

x w ^A^/ \ a • / / \ \ • / /^ / dre^- Xl > r Ai(r + S2 )Ai(r + S ' 2 ). (64) 

Ai(s 1 )Ai(s 2 )Ai(s' 1 )Ai(s 2 ) J 

We can check that Eq. (|64|) is invariant through uniform spatial translations. Next, from Pq 1 : q 2 (X\ , X 2 ) 
we can derive the distribution, Pq(X), of the relative Eulerian distance, X = X 2 — X\. It only depends 
on the relative Lagrangian distance, Q = Q 2 — Qi, through 

/oo 
dX 1 P Ql>Ql+Q (X 1 ,X 1 +X). (65) 
-OO 
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Fig. 9 (Color online) Left panel: The probability distribution Pq(X) of the Eulerian increment X, for three 
Lagrangian lengths, Q — 0.5, 1 and 2, from Eq. (|68jl . The integral over X > is smaller than unity as there 
is a n ad ditional Dirac contribution, Pg' ock 5(X), associated with shocks, with the weight p^ hock displayed in 
Fig. 1111 Right panel: Same as left panel but on a logarithmic scale. 



This gives 

" +io ° dsdsids 2 e ^Q-x)+(s 1 +s 2 )x/2+(s 1 -s 2 ) 2 /(iX) 



X>0: P Q {X) 



(27ri) 3 s 2 Ai(si)Ai(s 2 )Ai(si - s)Ai(s 2 - s) 
xl dre Xr Ai(r + s x )Ai(r + s 2 ), with 3?(s) < 0, (66) 



where the integration contour over s runs to the left of the pole at s = 0. Note that the expression 
(|66|) is similar to the result (|18|) obtained for the distribution of the Lagrangian increment Q over a 
fixed Eulerian interval X, except for the double derivative with respect to X and the factor 1/s 2 . This 
leads to the relationship between the distributions of the Eulerian and Lagrangian increments: 

d 2 d 2 

PQ(X) = — i Px{Q). (67) 



dQ 2 vv ' dX 2 

On the other hand, using the expression (| 1 1 1[) and the comparison with Eq. (|18p . taking the derivatives 
with respect to X in Eq. (|6"6")) gives the relationship 



P Q (X) = P x (Q) + 2^ J 



+loo dsdsids 2 e 'Vi-X) e -*.i.. a W 



(27ri) 3 s 2 Ai(si)Ai(s 2 )Ai(si - s)Ai(s 2 - s) 
x [-2sh SuS2 (X) - <P' SuS2 (X)h SuS2 (X) + h' Si . S2 (X)] . (68) 

We show in Fig. [9] the distribution Pq(X) obtained for three Lagrangian intervals Q. In a fashion 
similar to the Eulerian distribution Px(Q) shown in Fig. [TJ on large scales, Q ^> 1, the Lagrangian 
distribution Pq(X) is centered on Q, with cubic exponential tails on both sides as seen in (|69| below, 
whereas on small scales, Q <gC 1, it shows a monotonous decline. However, contrary to the Eulerian 
distribution, the Lagrangian distribution Pq(X) does not show an inverse square-root tail at low X 
as Pq(0) is finite. As can be seen in Fig. [9l the distribution Pq(X) given by Eqs.([66 | - ([68|) over X > 
is not normalized to unity as its weight decreases for smaller Q. Indeed, there is an additional Dirac 
contribution associated with shocks, of the form PQ lock S(X), where the Eulerian increment X is zero 
(all particles in the initial range [Qi,Q 2 ], of length Q — Q 2 — Qi, have merged into a single shock). 
Since the weight of this contribution grows at smaller Q, as can be checked in Fig. [11] below, the 
normalization of the regular contribution (|66[) decreases at smaller Q. 

On large scales, Q 1, from the expression (|6"6"|) or the relation (f6"T)) . we obtain 

Q > 1 : P Q (X) ~ T^X - Q), whence P Q (X) - P X (Q) for large Q and X, (69) 
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Fig. 10 (Color online) Left panel: The scali ng f unct ion Go(X) that describes the distribution of t he E ulerian 
increment X in the limit Q C 1, from Eqs. (|70[l - (|71[1 . The dashed line is the asymptotic behavior (|72[) . Right 
panel: Same as left panel but on a logarithmic scale. 



where the function Too was given in Eq. (|25p and shown in Fig. [3] Thus, at large scales the distribution 
Pq(X) is peaked around X — Q, with fluctuations of order unity that become increasingly small as 
compared with Q for Q — » oo. Moreover, it becomes identical to the Eulerian distribution in this limit. 

On small scales, Q C 1, using the expression (JTTTJ), we obtain a scaling similar to the one seen in 
Eg. ([50)) for the Eulerian distribution, 

Q<1: Pq{X)~QQq[X), withX>0, (70) 

and 

/+ioc i i _ 31 S2 (X) 

. Ioo »*iMiW K '' ( *- m -'«' ml ' (71) 

In particular, Eq. (|7ip gives the asymptotic behavior at large X: 

X»l: g (X)~ — ^ r^-mx-* 3 /^ (72) 
8Ai'(— u>iY 

We show the scaling function Go(X) in Fig.[Tnj We can see that it is finite at X = and is monotonically 
decreasing. 

The scaling of Eg. ([70)) is related to the fact that the system is described by a finite number of 
shocks per unit length, with masses of order unity (in terms of the dimensionless variables Q and X). 
Then, in the limit Q — > 0, the probability p^ hock that all particles in the interval [Qi,Q2], of length 
Q = Q2 — Qi, belong to the same shock goes to unity, and there remains a probability of order Q that 
the particles Q\ and Q2 belong to different shocks, in which case their Eulerian distance is of order 
unity (i.e. X ~ 1). This gives rise to the scaling (|70j) for this contribution associated with X > 0. 
This discussion shows that the scaling (|7D|) is less general than the scaling (|5U|) obtained for the small- 
scale Eulerian distribution, since it relies on the fact that shocks are well separated by distances of 
order unity. For instance, in the case of Brownian initial velocity, the scaling (|30[) is still satisfied but 
the property (|70|) is no longer valid. Indeed, in this case shocks are dense in Eulerian space and the 
typical Eulerian distance X between particles initially separated by the Lagrangian distance Q scales 
asX~ \fQ, as can be seen in [34[. 

As discussed above, in addition to the contribution associated with X > 0, there is a second 
contribution, of the form Pg hock 5(X), associated with the case where both particles Qi and Q2 belong 
to the same shock, whence X2 — X\. Its weight can be derived by computing the weight Pq{X > 0) of 
the contribution ([fjtj]) . Integrating Eg. ([fjtj|) . which gives two boundary terms at X = and X = +00, 
yields 

/W*>0)-I-i*- *(,)<„, ,73, 
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log(Q) 



Fig. 11 (Color online) Left panel: The probability p^ hock that a Lagrangian interval of size Q has collapsed 
into a single shock, from Eq. (|74|l . The dashed line is the asymptotic behavior (|75(l . Right panel: Same as left 
panel but on a logarithmic scale. 



where the integration contour over s runs to the left of the pole at s = 0. The comparison with 
Eqs.([30 l) -([3l ]) gives the relation 

J2 poo 

pshock = MQ) ^ whence pshock = / dQ , {QI _ Q) MQl y (?4) 

a V ' JQ 

In particular, using the results of section Eq. (fT4")) gives at once the asymptotic behaviors 

Q^O: P^ hock - 1-JF (0)Q~ 1-1.674Q, Q -» oo : P^ hock - 32^Q- 3 / 2 e^ lQ - Q3 / 12 . (75) 

Thus, we recover the fact that the probability, Pg hock , for two particles of initial Lagrangian distance 
Q, to belong to the same shock, goes to unity for Q — > 0, with a linear deviation so that Pq{X > 0) ~ 
Fo(Q)Q, in agreement with the scaling ([70f . 

We display in Fig. [TT] this probability Pg hock , which clearly shows its steep falloff at large Q. Again, 
the cubic exponential decay can be understood from the same arguments, as those used for the tails 
(fTTjl or (I109P of Eulerian distributions. Note that the formation of a single shock of strength q = q2 — qi 
would be associated with a velocity difference vi — v\ = —q/t, rather than with the mean velocity vq 
over the length q as was the case for these Eulerian distributions. However, this difference does not 
give well-defined results for the initial white-noise velocity field. Again, this is due to the fact that the 
dynamics is governed by non-local processes, that is, one cannot obtain behaviors such as ([75*)) from a 
local analysis (i.e. a local Taylor expansion) of the initial velocity field. This clearly follows from the fact 
that, at any time t > 0, matter has gathered in a series of discrete shocks, which has strongly modified 
the velocity field: the latter has been regularized by the balance between the (infinite) different sign 
velocities of neighboring particles over lengths of order (2Dt 2 ) 1 ^ 3 and become strongly non-Gaussian, 
as seen in sections lB.ll and 14.21 However, to recover the tail (|75|) . we can split the interval q into two 
equal parts, and note that at least one of the two mean initial velocities v\ and V2 of both intervals 
must be of order q/t, which leads back to the cubic exponential tail ([75]) . 



6.3 Mass function of shocks 

We briefly note here that the mass function of shocks, n(m), can be derived from the shock probability 
pshock g^vrdied in the previous section. Here we define n(m)dm as the mean number of shocks, per 
unit Eulerian or Lagrangian length (both functions are identical since on large scales X = Q up to 
fluctuations of order unity, as seen in sections 14. 21 and 16 . 2p . with a mass in the range [m,m + dm\. As in 
section [5j since we consider a uniform initial density po, we have m = paq for the mass associated with 
a Lagrangian interval q. Then, the probability, Pq hock , that two particles of initial Lagrangian distance 
q belong to the same shock, can be obtained by counting the number of shocks of mass m > poq, each 
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shock giving rise to a factor (m/po — q) as q\ may be located within the distance (m/po — q) from its 
left boundary. In terms of dimensionless variables this reads as 

pshock = / dM (M — Q)N(M), with n(m) = =N(M) and M = — . (76) 

Jq * PoT Pol 

Then, using Eq. (|7"4")) we obtain at once 

= MM). (77) 

Q=M 

From the expression (j32|) we can check that we recover the result of [l3| , who directly derived the shock 
mass function from the geometrical construction ^ without considering the Lagrangian displacement 
field. This provides a useful check of the computations performed in section [6~^1 within a Lagrangian 
framework. The asymptotic properties of the shock mass function can also be read from Eq. (|55|) . see 
also [!, H, HH . Moreover, the integral properties 

/>oo />OC 

/ dMN(M) = F a (0) ^ 1-674, / dM M N(M) = -^(0) = 1, (78) 
Jo Jo 

ensure that mass is conserved and that there is a finite mean number of shocks per unit length 
(~ 1.674 shocks in the mean, in units of X and Q). Many more properties of shocks, such as their 
n-point multiplicity functions as a function of their mass and velocity, can be found in [l3j ] . 



^2^>shock 

N(M) = 



7 Small-scale heuristic approach for general initial conditions 

The results obtained in the previous sections were derived from exact computations, based on Eqs. (|12p 
and (|13[) . In this section, using a heuristic approach that assumes that small-scale properties are 
governed by shocks, or point-masses in higher dimensions, we discuss how small-scale scalings obtained 
for 1-D white- noise initial conditions would extend to generic initial conditions and higher dimensions. 



7.1 General 1-D case for the distributions of Lagrangian and velocity increments on small scales 

As explained in section [431 the scalings Px(Q) ~ XJ 7 (Q) and (Q v ) cx X are due to the presence of 
shocks and as such they apply to a large class of initial conditions [l||, [3(| • Then, the distributions of 
the Lagrangian increment q and of the velocity increment v over the distance x are determined by the 
one-point distribution of shock strength [30] , and they factorize as 

x^O, q > 0, v<x : p x (q)~xn(q), p x (v) ~ xt n(x — v) ~ xtn{— v), (79) 

where n(q) is the mass function of shocks, that is n{q)dq is the number of shocks of strength q per 
unit length. Here we used the dimensional variables x and q because the power of time that appears in 
the relevant scaling variables X and Q depends on the initial conditions. For instance, for a power-law 
initial energy spectrum, Eo(k) cx k n , with —3 < n < 1, we would have X cx x/i 2 /( 3+n > 0. For 
initial energy spectra that are not a power law there may not exist scaling variables such as (|5|) (for 
instance, for a smooth spectrum one may expect a time-dependent effective exponent n(t)) but as soon 
as shocks are present one still has scalings of the form {q v ) ~ x at small distance for v > 1. Note that 
for v = 1 we always have the exact relation (q) = x, because of the conservation of matter, but the 
linear scaling over x does not always extend to powers < v < 1, as for the white- noise case studied 
in the present paper, see Eq. (|34[) . For instance, for Brownian initial velocity it only extends down to 
v = 1/2 [H, [13]. Indeed, as v decreases the moment {q v ) becomes increasingly sensitive to low-density 
regions, characterized by low Lagrangian increment q, so that shocks are no longer dominant (for the 
white-noise case they remain dominant down to v = + because shocks are separated by voids where 
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the Lagrangian increment q is exactly zero). It can be useful to introduce the moment generating 
function W x (s), defined by 

00 (—a\n poo 

^ = E -^~ {qn) = J dq ^ sq ~ ^ Px[ql m 

whence 

r +io ° ds 

Px(g)= / ^e'« for q>0, (81) 

where we assumed that all moments are finite and uniquely determine the function \P x (s) (note that 
adding a constant to ^(s), so that ^(0) 7^ 0, does not contribute to p x (q) for q > as it only yields 
a Dirac S(q)). From the small-scale scaling ([751) we obtain 

x^O: W x {s) ~ xW(s) and (<f) - x (-1)™ <F (n) (0), with W(s) = f dq (e^ 39 - l) 71(g). (82) 

./o 

For the case of white-noise initial velocity studied in this article, using the dimensionless variables X 
and Q, we have N(Q) = F Q (Q) and W(s) = F (s) - F Q (O^as shown by Eqs. (J3DI»-(EI1» and Eq.^TTJ). For 
the case of Brownian initial velocity, using results from [34| , we have 

N(Q) = —= Q~ 3/2 e~ Q , W{s) = 2(1 - y/i + s), for Brownian initial velocity, (83) 

where we used the relevant scaling variables of the form Q oc q/t 2 (here Eo(k) oc fc~ 2 ). We can check in 
[34| that N(Q), given in Eq. (|83|) as the shock mass function, also describes the probability distribution 
Px(Q) through Px{Q) ~ XN(Q) on small scales, as explained above in ([75|) . Note that this remains 
valid even though shocks are no longer isolated but dense in Eulerian space. Indeed, if we select shocks 
above a small finite mass threshold m*, the latter are again isolated so that the previous arguments 
apply, and smaller shocks only modify the low-q tail of the distribution. Indeed, the previous arguments 
hold for the limit x — > at fixed q, or more precisely above a cutoff g_ (x) that goes to zero faster 
than x, so that the scalings (q u ) oc x for v > 1 do not depend on the behavior of the distribution 
Px(q) over this low-g domain. Thus, the functions n(q) and ^(s) do not necessarily describe the actual 
distribution p x (q) down to q = for a finite x. As seen above, for the case of white-noise initial velocity, 
the scaling function !Fo(Q) of Eqs. (|3"0|) - (|3"Tjl actually applies down to Q = + , as it only misses the 
Dirac contribution (I108j) (and the very high-Q tail (|23|) which is repelled to infinity). However, for the 
case of Brownian initial velocity for instance, the scaling functions (f83|) only apply to Q > X 2 and 
they miss a \ow-Q cutoff of the form e^ x2 ^ [34]. 



7.2 General 1-D case for the distribution of overdensities on small scales 

As seen in sect. l5.ll since the overdensity 77 is also given by the ratio 77 = q/x, its probability distribution 
is related to the distribution of the Lagrangian increment q through p x (v) — x Px{q)- Then, from 
Eqs.fTJ)-® we have 

r +ico ds — 

x^O, i]>0: p x (n) ~ x 2 n{x7 1 ) = x 2 / — e sx,1 V(s). (84) 



2?ri 

If we define the cumulant ratios S n and their generating function <p x (y) as in Eqs.([48l)- (T49|) . 

Sn(x) = f^ and v x {y) = Y J {-l) n - 1 S n (x) y -, (85) 

W )c n=1 

using the property that in the small-scale limit, x — > 0, we have for 77 > 1 the asymptotic relationships 
(V n )c ~ (q n )/% n > we obtain 

x^O: <p x (y)~m with (p{y) = -W"{Q)W\ ) , (86) 

\# (0)7 
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and _ 

x^Q: Px ( v ) ^ - J +l °° ^-^- e y^^^(y) with (if) c = *M. (87) 

We can check that this agrees with Eq. (|151) for the white- noise initial velocity studied in this article. 
For the case of Brownian initial velocity [34(, we have, in agreement with (|83[) . 

X — > : Px(v) ~ \ j — r/~ 3 / 2 e~ Xri and (p{y) = y/T+~2y — 1, for Brownian initial velocity, 

V 7T 

(88) 

where X oc x/t 2 is again the scaling variable relevant to that initial condition. For the Brownian case 
it happens that the expression ([85]) for the density cumulant generating function is actually exact for 
all scales, but in the general case, as discussed below Eq. (|83]) , it only applies to the small scale limit at 
hxed density contrast. More precisely, although Tq describes the overdensity probability distribution 
down to r\ — > + through Eq. iP^l) in the case of white-noise initial velocity, in the general case it only 
applies above a cutoff ry_ (a;) that decreases with x (for instance for the Brownian initial velocity we 
have r/_ (x) oc a;). 



7.3 Multifractal formalism in D dimensions 

These properties can be extended to higher dimensions D through a heuristic multifractal formalism 
0, QHH, without going through the inverse Lagrangian map x^q. Thus, let us assume that the 
overdensity within a spherical cell of radius I centered on x scales for I — > as %(x) ~ i a for points 
x e V a C R" , with dimD a = F(a). Then, we may write 

1^0, v>Q: (r#) = J™ 6tiTfp t (Ti)~ Jdaf a+D - F ^p*(a), (89) 

where p*{cx) gives the weight of the various multifractal exponents, and we used the fact that the 
probability for a sphere of radius t to encounter an object of dimension F scales as £ D ~ F for I — > 0. 
Using a steepest-descent argument, we obtain the small-scale exponents 7„, 

£->0: {rft) ~ t" 1 " with j v = - mm[va + D - F(a)} = max[F(a) - va - D]. (90) 

a a 

Thus, the exponents 7„ and F(a) are related by a Legendre transform. The fractal scaling exponent a v 
that is associated to 7„ through ([90|) is the abscissa of the first-contact point of the curve F(a) with the 
family of straight lines, va + c, moving downward from c = +oo, in a fashion similar to the geometrical 
construction associated with ([3]). In particular, the exponents 7„ only probe the concave hull of F(a) 
[Til ]. Since the matter density is positive the scaling exponents a„ are restricted to a > —D. This 
lower bound corresponds to Dirac density peaks (i.e. massive points), which have a zero dimension, 
F(—D) = 0. On the other hand, the constraint associated with the conservation of matter, (rfi) = 1 
whence 71 = 0, ensures that the curve F (a) is below the straight line a + D, which runs through the 
point {— D, 0}, and has at least one contact point with this line. Then, we can see that, as soon as 
isolated Dirac density peaks have formed, with a finite probability per unit volume, the first-contact 
point between F(a) and the family of straight lines va + c with v > 1 is the point {—D, 0} (for v = 1 
there can be other additional contact points), which gives j u = (y — \)D for v > 1. For instance, in 
three dimensions D = 3, we generically expect to first form "Zeldovich pancakes" [23], that is sheets 
with a finite surface density, that intersect to form filaments (i.e. lines of finite line density), which join 
to form point-like masses (nodes). This corresponds to objects of fractal exponents and dimensions 
{ — 1, 2}, {—2, 1}, and {—3, 0}, along the line {a, F = a + D}, that all contribute to 71 = while only 
point-like masses contribute to 7„ = "&(y — 1) for v > 1. Then, we obtain 

n > 1 : (r) n ) - r {n - 1)D whence (rf ) c - r^~^ D , S n (£) ~ t Q and tp t (y) - <p{y) for t -»• 0, 

( 91 ) 

where the limiting generating function ip(y), reached in the small-scale limit, no longer depends on 
£. Thus, for stochastic initial conditions, where we generically expect the formation of isolated Dirac 
density peaks in finite numbers per unit volume, we obtain the scaling (|9ip and the first expression in 
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(f87| for the density probability distribution, where in dimension D we have (t] 2 ) c ~ £~ D . As explained 
above, within this heuristic multifractal formalism the property (f91"j) holds independently of the form 
of the fractal spectrum F(a) over a > —D, as soon as there is a finite density of point-like masses 
which gives rise to the fractal exponent {a = —D, F = 0}. The computation of the complete spectrum 
of fractal dimensions, that is necessary for the study of exponents v < 1 in Eq. |89|) . is a difficult task 
and in the general case there can be a continuous rather than discrete spectrum (as discussed below 
in a phenomenological fashion for the case of Brownian initial velocity, where shocks are dense) . 

Again, the scaling function (p(y) and the distribution it defines through (|87|) only apply above 
a density threshold rj-(£) that may only show a very slow decrease with £. Within the multifractal 
formalism ([55)1 the behavior at small densities, below this threshold ry_ (£) , depends on the properties of 
the curve F(a) to the right of the point {— -D,0}, see [3, Hl|. For the one-dimensional case with white- 
noise initial velocity studied in this article, since shocks form isolated density peaks amid empty space, 
the curve F(a) is reduced to the single point {—1, 0} so that the scalings (rf) ~ X x ~ v apply to all v > 0, 
in agreement with the second relation in (|46[> and the fact that the distribution Px(v) ~ X 2 Fo(Xr)) 
in (pro]) applies downto rj = + , as seen in the previous sections. 

For the one-dimensional case with Brownian initial velocity, we have the bifractality (q u ) ~ x for 
v > 1/2 and {q") ~ x 2u for v < 1/2 PHI. This leads to {rf) ~ x 1 ^ for v > 1/2 and (17") - x u for 
v < 1/2. In terms of the multifractal formalism, this would be interpreted as a second point {1,1}. 
This may be understood as follows. For these initial conditions, the shock mass function diverges at 
small masses as n(m) oc m -3 / 2 and shocks are uncorrelated and dense in Eulerian space [28l. [29l la. l34j . 
Then, choosing a small finite mass threshold m*, the set of shocks of mass larger than m* gives a 
population of isolated point-like masses that gives rise to the fractal exponent {—1,0}. On the other 
hand, if we choose a random Eulerian interval of size £, it contains in the mean m~ 1 ^ 2 £ shocks of mass 
in the range [m, 2m ]. Taking m ~ £$ we obtain that i) for any /3 > 2 an interval of size i contains of 
the order of f-^P/ 2 shocks of mass in [m, 2m], which leads to an overdensity larger than r\g ~ £^/ 2 , and 
ii) for any [3 < 2 an interval of size £ contains with a probability of order £ X ~P/ 2 at least one shock 
of mass in the range [m, 2m], which leads to an overdensity larger than r\i ~ £^~ x . The first point 
i) leads to the multifractal point {1, 1}, and the second point ii) leads to the points {(3 — l,/3/2} for 
< f3 < 2, that is to the segment joining the points { — 1, 0} and {1, 1}. Therefore, we obtain in this 
case the multifractal spectrum F(a) = a/2 + 1/2 with — 1 < a < 1, see also [l[ for more rigorous 
discussions. More generally, in one dimension the scalings obtained for v < 1 and the low-density tail 
are related to the low-mass tail of the shock mass function. 

We can note that the small-scale limit of finite ratios S n and generating function ip x (y), as in 
(|86]l . also corresponds to the "stable-clustering ansatz" introduced in the cosmological context as a 
phenomenological model for the highly nonlinear regime (27j . There, it was derived by assuming that 
on small physical scales, after nonlinear collapse and gravitational relaxation, overdensities decouple 
from the Hubble expansion and keep a constant physical size [13]. It can also be associated with a 
multifractal formalism, where the moments of the density with v > 1 would be governed by a single 
fractal exponent as in (|91|) . which however would not be associated with point-like masses but with 
structures of exponent a ~ 1.3 and dimension F — 3 — a ~ 1.2 [l,[3l[. However, contrary to the Burgers 
dynamics, this behavior may not be exactly reached on small scales for the gravitational dynamics, 
as the coefficients S n still appear to show a weak dependence with scale On the other hand, the 
Burgers dynamics itself is also known as the "adhesion model" in this cosmological context fl8l. [3al , 
where it provides a good description of the large-scale filamentary structure of the cosmic web [23| • It 
is not clear whether the reasonably good match of the "stable-clustering ansatz" could be understood 
from the the exact scaling (|9Tj) achieved in the small-scale limit by the "adhesion model", since the 
nonlinear structures are different (point-like masses as opposed to extended halos) and no detailed 
comparisons have been performed yet in terms of the ratios S n themselves. 



8 Conclusion 

We have obtained in this article some equal-time properties of the Burgers dynamics, in the inviscid 
limit for white-noise initial velocity. In agreement with previous works, the initially singular distri- 
butions are regularized as soon as t > 0, but little power is transfered to large scales. Thus, the 
distributions of the fluctuations of the Lagrangian increment, q, and of the velocity increment, v, 
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around their means, have a finite limit in the large-scale limit x — > 0. We recover the characteristic 
cubic exponential tails associated with white- noise initial conditions. Voids lead to an additional Dirac- 
type contribution to these distributions, that also decays as a cubic exponential at large scales and is 
preceded by an inverse square root tail with a weight of the same order. On small scales, where the 
probability to be within a void goes to unity, the regular part factorizes as X Tq{Q), which corresponds 
to the probability to contain one shock of strength Q. In particular, the scaling function Tq{Q) is also 
the mass function of shocks. This leads to the standard linear scaling with x of the velocity structure 
functions at small scale, due to shocks. 

Next, we have derived the distribution of the density within intervals of size x. It presents simi- 
lar properties to those obtained for the Lagrangian increment, and exhibits the corresponding large- 
scale and small-scale scalings. In particular, at small scales this gives rise to the scaling hierarchy for 
the density cumulants known as the "stable-clustering ansatz" in cosmology. Here it is due to the 
presence of shocks. We also obtain the density two-point correlation and power spectrum, with the 
high-wavenumber constant asymptote associated with shocks. 

Turning to the Lagrangian displacement field, associated with a description of the dynamics in 
terms of Lagrangian coordinates, we have obtained the distribution of the Eulerian increment x for a 
given mass poq. On large scales the Lagrangian distribution p q (x) becomes identical, at leading order, 
to the Eulerian distribution p x (q). On small scales there is also a factorization of the form QQq(X), 
but this is less general than for the small-scale Eulerian distribution since it only applies to initial 
conditions such that shocks are isolated, that is initial energy spectra with — 1 < n < 1, whereas 
the Eulerian factorization remains valid for the whole range —3 < n < 1. Contrary to the Eulerian 
distribution, the Lagrangian distribution p q (x) does not show divergent tails as it remains finite for 
x — > 0, but there is again an additional Dirac contribution, which is now due to shocks. 

Finally, within a heuristic approach we have discussed how these small-scale properties generalize 
to other initial conditions and give rise to a universal scaling for the distribution of the Lagrangian 
increment (and of the velocity increment) above a lower cutoff q~(x), that goes to zero faster than x 
in a fashion that depends on the initial conditions. A heuristic multifractal formalism allows to extend 
these results to higher dimensions. It generically leads to a universal scaling hierarchy for the density 
cumulants in the small-scale limit, that is governed by point-like masses. This also corresponds to 
the "stable-clustering ansatz" introduced in the cosmological context, where the Burgers dynamics is 
known as the "adhesion model" and is used to describe the large-scale cosmic web. 

The results obtained in this article may prove useful to test approximation schemes devised to 
handle other initial conditions or closely related dynamics, such as Navier-Stokes turbulence or grav- 
itational dynamics, where no exact results are available, as in [l]|, H3, HH . In this respect, the case 
of white- noise initial velocity studied here would present a severe test for non-perturbative methods. 
Indeed, the initial energy spectrum is so "blue" that nonlinear effects are dominant up to the largest 
scales, x — > oo, and perturbative expansions already encounter ultraviolet divergences at leading or- 
ders. This implies that alternative approaches must be able to take into account shocks, as for the 
steepest-descent methods presented in [331 ] . Another interesting feature of the case of white-noise ini- 
tial velocity studied in this article is that it shows a density power spectrum that displays two different 
large-scale and small-scale regimes, as for the gravitational dynamics in the cosmological context, but 
can still be computed exactly. 



A Transition kernel with parabolic absorbing barrier 

For the white- noise initial conditions ((5}, the process q t—> tpo is Markovian and a key quantity is the conditional 
probability density K x , c (qi, q2, 1P2) for the Markov process ipo(q), starting from at 31, to end at ip2 at 
32 > 31, while staying above the parabolic barrier, tpo(q) > Vx,c(q), for qi < q < 32. We briefly recall here its 
derivation, obtained in [TST ]. using our notations. It obeys the diffusion equation 

d D d 2 

<72 > qi ■■ Q^K x , c {qi,ipi;qa,i>a) = — -^K x ,c{qx,i>v,q2,i>2) (92) 

over the domain ip > V x , c (q), with the initial condition at 32 = Si, ^Cb,c(si, ipi', Si, V'2) = S(tp2 ~ V>i), and the 
boundary conditions, K x<c (qi, tpi; 52, 1P2) = for ip\ = Vx,c(qi) or 1^2 = V x , c {q2)- The kernel associated with 
the propagation towards the left side, 52 < Si, is obtained from the parity symmetry 



52 < Si : K x , c (—q 1 ,ipi;—q2,ip2) = K- X<c {qx, ipV, 32, V^)- 



(93) 
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In terms of the dimensionless coordinates (JSj the kernel K x ,c can be written as 

K x , c (Qi^i;Q2,^)di, 2 = eU*-**- W 1 -^)^-(Q 2 -^ 3 /3+(Q 1 -x)3/3 G(r . ri) ra) dr2> (94) 



where we defined 



T — Q 2 ~ Qi, n = 2 



ft + W«-*> a -c 



(95) 



and the propagator G obeys the Schrodinger-like equation 

dG _ d 2 G 

-3- = -r 2 G + — - 
dr an 



— = ~r 2 G + 10 - over r > 0, r > 0, (96) 



with the initial condition G(0; 7-1,7-2) = S(r 2 — r\) and the boundary conditions G(T;ri,r 2 ) = for n = or 
72 = 0. This reduced propagator G can be solved as [T6l. [l3j] 

G(r; rr, r 2 ) = g (97) 

where — u>h are the zeros of the Airy function Ai(:r) (in particular, wi ~ 2. 338) . Thus, G(r;ri,r 2 ) is symmetric 
over {ri,72}, and it also obeys the backward equation (compare with Eq. (|96l) ) 



dG „ d 2 G 

or orf 



= —ri G + -t^j- over r > 0, r > 0. (98) 



Next, it is convenient to introduce the probability density, E XtC (qi, ipi; q 2 , ip 2 ; q)dqdcdifi 2 , for the curve ij>a(q), 
starting from t/>i at 51, to end at ip 2 at q 2 > gi, while staying above the parabolic barrier P IiC , and with a last 
excursion below V x ,c+dc in the range [q, q + dq]. From the definition of the kernel K x . c , it reads as 

E x , c (qi,ipi;q 2 ,ip 2 ;q) — lim / [K XlC (qi, ipi; q, ip) — K X)C+ s c (qi, ipv, q, tp)] K x>c (q, V"! 32, ipn). (99) 



dq Sc^o 5c 

Using a Taylor expansion and integrations by parts, this yields [l3| 

E XlC (qi,ipi;q 2 ,ip 2 ;q) = ^ 9 ^*' c (qi,ipv,q,ip) d ^' c (q, %l>\ q 2 , fa] 
2 oip 2 oipi 



(100) 



which gives in terms of the reduced propagator G introduced in Eq. (l94|) 

E x Aqi,^,q2,^;q) = ^ e (Q 2 -x),,-(Q 1 -x), 1 -(Q 2 -x)3 /3+ (Q 1 -x)3 / 3 

x ^(g-Qrsn, 0)^ r (Q 2 -Q; 0,r 2 ). (101) 

As could be expected, the expression (|10ip shows that the probability density E x , c depends on the behavior 
of the propagator G close to the boundary r = at one end. This corresponds to the contact point of abscissa 
q between the curve tpo and the parabola V XyC that is involved in the definition of E x . c . 

B Eulerian distributions 

B.l One-point distributions p x (q) and p x (v) 

Substituting Eqs. (|101[l and (I97|l into Eq, (|10[l . one obtains in terms of the dimensionless variables (|8jl. 

Px(Q) = J(X - Q) J(Q - X) and P(V) = J(V)J(-V), (102) 
where we used the relation X = Q + V , with [la ] 



j(u) = lim r dre -(-") 3 /3+(-^y e -^M!^ = r°°ii^ 

— Jo Ai'(-w fc ) 7_ loo 2^iAi 



A' = l 

From the asymptotic behaviors of the function ^(u). 



(103) 



it — > +00 : J" (it) ~ 7-777 T> anc ^ ^ or u — * — 00 ' ^( u ) ~ — 27te" ,/3 , (104) 

Ai ( — u)i) 

one obtains the asymptotic behavior of the distribution of the velocity V (and of the Lagrangian coordinate 
Q = X-V)^. 
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0.5 1 1.5 

X = x/(2DL 2 )>- /3 



-1.5 -1 -0.5 0.5 

log(X) = log[x/(2DL 2 ) 1 - /a ] 



Fig. 12 (Color online) Left panel: The probability P% that an Eule rian interval of size X is empty, that is, 
that the Lagra ngian increment Q over this interval is zero, from Eq. (|108[) . The dashed line is the asymptotic 
behavior (|109|l . Right panel: Same as left panel but on a logarithmic scale. 



B.2 Two-point distributions p Xl ,x 2 (<Zi) 92) p XltX2 {v\ 1 V2) 

We first consider the case i) of section I3T21 when the two first-contact parabolas V X1 , C1 and V X2 ,c 2 have two 
different contact points gi and qi with the curve ipoig) (and there is at least one shock in the interval [xi, X2] 
since the map x 1— > q is constant outside of sho cks [2a. Then, noting the abscissa of the intersection 

between both parabolas, in a fashion similar to (I10|l we can write this contribution to p x±:X2 (qi, 52) as 

(105) 



P?i,z 2 (91,92) = lim / dcidczdip*dxl>+E XUCl (q-,a;q*,ip*;qi)E X2tC2 (q* ) xl>y,q + ,ip + ;q2). 



Substituting the expression (|10ip of the kernel E X)C gives, in agreement with [l3l ] . the expression (|12|l . where 
we introduced the function Iri defined by 



Hx 1 ,x 2 (Qi,Q 2 ) =2{X 2 -Xx) / dr, / dQ. e { 

JO J Qi 

x-5— (<?* -Qi;0,r») ^— (Q2 - Q,;r,,0). 



9ri 



9r2 



(106) 



We can check that the function Hx 1 ,x 2 (Ql, Q2), whence the distribution Pj^ x 2 ( ( ?i> < 2 2 )> are invariant with 
respect to uniform translations of Xi and Qi, in a gree ment with the statistical homogeneity of the system. 

We now consider the second case ii) of section [3721 when the two parabolas intersect at the common point 
qi =qi, and we can write this contribution to p Xl: 

p7 1 ,x 2 (qi,Q2) = 8(q 2 - qi) lim lim / dcidi/>*d?/>+ E Xl)Cl (<?_, 0; q,,i^,;q\) K X2 , C2 (q,,^,\q + ,%f} + ). (107) 
This gives in terms of dimensionless variables [TTJ] the expression (|13p . 



B.3 Probability Pji of empty intervals of size X 

Using Eq, (|103|l . we can write the second term in (|17p as [L| 

Px{Q) = 5{Q)P°x with Pl = Jfr^/ ^ T77TA7T71 ' 

VA J_ ioc (2tti) 2 Ai(si)Ai(s 2 ) 

This yields for the probability P% to have a vanishing Lagrangian increment the asymptotic behaviors 



X -> : P¥ 



and for X — > 00 : P x 



Ai'(-Wi) 



£ X" 1/2 e -"i^-* 3 /i 2 



(108) 



(109) 



Since Eulerian intervals with Q = have a zero matter conte nt, P ° is also the probability for an interval of 
size X to be empty, in agreement with the result obtained in [l3l ] for this void probability. We compare this 
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probability P% with its asymptotic behavior (| 109[) in Fig. 1121 see also [T3| . The cubic exponential tail (| 109[) may 
be understood using the same arguments as those used for the tail of the velocity distribution |TT} discussed 
above. Thus, for the Eulerian interval of size x to be empty, its initial matter content must have traveled by a 
distance of order x, which requires a mean velocity over this interval of order v ~ x/t. Again, since the initial 
Gaussian velocity over scale x is vo(x) = (ip2 -~ipi)/x, with a variance a% () (x) = D/x, this yields the probability 
^ e -Wt) 2 /<^ M ^ e -^ 3 /(ot 2 ) ; which gi ves b ac k tne cu bic exponential tail JT09J. 



C Laplace transform of the product of two Airy functions 

We recall here the results obtained by pl|] for the integral over two Airy functions that appears in Eq. (|18|) . 
Thus, if we define /(r) and g(x) by 

/(r)=Ai(r + ai)Ai(r + a a ), g(x) = / dre"7(r), (110) 

Jo 

the Laplace transform g(x) can be integrated as 

g(x) = -±=e*™M-e*™M dye-*^^h sl , S2 (y), (111) 

Jx 

with 
and 

ft (,)_/("), AO) , r(0)-2( Sl+S2 )/(0) /( 3 )(0)-2( Sl+S2 ),f(0)-2/(0) 
S1,s 4 4 4x 4a; 2 

The first term in the right hand side of Eq.{TTTJ gives the large- a; behavior of g(x), up to terms of relative 
order e~ x /12 . 
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